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Abstract: We study the mismatch between a full calculation of non-global single- 
logarithms in the large- N c limit and an approximation based on free azimuthal averaging, 
and the consequent angular-ordered pattern of soft gluon radiation in QCD. We compare 
the results obtained in either case to those obtained from the parton showers in the Monte 
Carlo event generators HERWIG and PYTHIA, with the aim of assessing the accuracy of 
the parton showers with regard to such observables where angular ordering is merely an 
approximation even at leading-logarithmic accuracy and which are commonly employed for 
the tuning of event generators to data. 
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1. Introduction 

An important class of theoretical predictions in QCD fall under the banner of "all-order" 
calculations. This refers specifically to predictions for those observables that receive loga- 
rithmic enhancements at each order of perturbation theory, which threaten the convergence 
of the perturbation expansion in important regions of phase space. A classic example is 
event-shape distributions, where studying an observable close to its Born value (such as 
the distribution of the thrust variable 1 — T near T = 1) results in generating terms as 
singular as cXsjhf hi 2n-1 (l — T) in the perturbative prediction, which can render all orders 
in a s equally significant The origin of these logarithmic enhancements is the singu- 
lar behaviour of the QCD emission probabilities and their virtual counterparts in the soft 
and/or collinear kinematical regions. These singularities coupled with the nature of the 
observable (where measuring close to the Born value constrains real emission but not the 
purely virtual terms) lead to the appearance of large uncancelled logarithmic contributions 
in the fixed-order perturbative results. 

There exist two main approaches to deal with such logarithmic enhancements at all 
orders. The first is the method of analytical resummation where insight on the QCD 
multiple soft-collinear emission probabilities and analytical manipulations of the phase 
space constraints are carried out 1 so as to obtain a result that resums the large logarithms 
(for those variables that satisfy certain conditions ensuring they can in fact be resummed 
@) into a function which can be expressed in the form 



There exist a variety of formal approaches designed to achieve these goals all of which embody the 
physics that we outline here. 



£(F) = exp [Lgi(a s L) + g 2 {a s L) + a s g 3 (a s L) -\ ] 
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where Lgi, g2, etc. are functions that are computed analytically 



2 and V is a generic 



event shape, e.g. 1 — T. The function Lgi, if non-zero, represents the leading or double- 
logarithmic contribution (LL), since it contains an extra power of L relative to the power 
of a s , i.e. C(a™L n+1 ). g^ is the single-logarithmic or next-to-leading logarithmic (NLL) 
contribution containing a logarithm L for each power of a s , (D(a™L n ), etc. We especially 
note that if the function gi is zero (as in the case of the interjet energy flow observable 
we shall study in detail here), the single-logarithmic function g2 contains the leading loga- 
rithms. The function a s g^ contains an extra power of a s relative to the power of L and is 
next-to-next-to leading logarithmic (NNLL) if gi is present and next-to-leading logarith- 
mic (NLL) otherwise. In the limit V ^ 0, S(V) has a physical behaviour as opposed to its 
expansion to any fixed order, which is divergent as we mentioned. This expression, which 
is valid at small V, can then be matched to exact fixed-order estimates that account for 
the large-1/ region, so as to give the best possible description over the entire range of V. 

Another possible approach to studying such observables is provided by Monte Carlo 
event generators amongst which the most commonly employed are HERWIG [Q, |5| and 
PYTHIA |], [fl, g|. We note that these programs are of far greater general utility than 
the study of the observables we will discuss here, providing simulations of complete QCD 
events at hadron level and representing perhaps the most significant physics tools in current 
high-energy phenomenology. The parton showers contained in these event generators aim 
to capture at least the leading infrared and collinear singularities involved in the branching 
of partons, to all orders in the large-Ac limit. One may thus expect that the dynamics 
that is represented by the parton shower ought to be similar to that which is used as ana- 
lytical input in QCD resummations at least on the level of the leading (double) logarithms 
involved. 

For several observables a correspondence between the Monte Carlo parton shower 
and the matrix elements used in analytical resummations is in fact clear. Considering, 
for example, final-state radiation, parton showers evolve due to parton emission with the 
branching probability V satisfying || [H], |ll| 



where te^ax * s the maximum k 2 accessible to the branching and is a cut-off regular- 
ising soft and collinear singularities. The above result, with P(z) being the appropriate 
Altarelli-Parisi splitting function relevant to the branching, captures the soft (z — * 0) and 
collinear (k 2 — ► 0) singularities of the emission. Virtual corrections (and hence unitarity) 
are incorporated via the Sudakov form factors A(/c 2 , k^). 

An essentially similar form is employed for the purposes of most analytical resum- 
mations where the probability of emitting several soft gluons is treated as independent 
emission of the gluons by the hard partons which for simplicity, in the rest of this paper, 
we take to be a qq pair. The probability for emitting a soft and/or collinear gluon is the 

2 We include in the category of "analytical" the semi-analytical approach of Ref. Q where analytical 
observations are exploited such that Q2 can be calculated numerically in an automated fashion for several 
observables. 




dV 



^(^maxi ko) 



(1.2) 
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very form mentioned above and the virtual corrections are included as in the Sudakov 
factor. This independent-emission or probabilistic pattern (which stems from the classical 
nature of soft radiation) suffices up to next-to-leading or single logarithmic accuracy for a 
large number of observables. Thus it is natural to expect that at least as far as the double- 
logarithmic function g\ is concerned, it would be accurately contained within the parton 
shower approach, although it cannot be separated cleanly from the single-logarithmic and 
subleading effects generated by the shower. Beyond the double logarithmic level one ex- 
pects at least a partial overlap between the parton shower and the analytical resummations, 
where the degree of overlap may vary from observable to observable and depend on which 
hard process one chooses to address. The state of the art of most analytical resummations 
is next-to-leading logarithmic, i.e. computing the full answer up to the function gi- Monte 
Carlo algorithms such as HERWIG are certainly correct up to g\ and perhaps in certain 
cases gi accuracy (while being limited to the large- iV c approximation) but not beyond (see, 
e.g., the discussion in |l2| ). 

As we mentioned, the event generator results do not explicitly separate leading loga- 
rithmic from next-to-leading logarithmic or subleading effects (e.g. those that give rise to 
53 and beyond) and, moreover, parton-level Monte Carlo results include non-perturbative 
effects that arise via the use of a shower cut-off scale, i.e. k$ in Eq. (|l.2| ). Erom the 
point of view of having a clean prediction valid to NLL accuracy that can be matched to 
fixed-order and supplemented by, for instance, analytically estimated power corrections, 
one would clearly prefer a resummed calculation. This is not a surprise since these calcula- 
tions were developed keeping specific observables in mind unlike the event generators which 
have a much broader sweep and aim. It is thus not our aim to probe event generators as 
resummation tools in themselves but rather to consider the logarithmic accuracy to which 
perturbative radiation may be generically described by a parton shower of the kind to be 
found in HERWIG or PYTHIA, for different observables. 

The above is particularly important since it has been pointed out relatively recently 
that for a large number of commonly studied observables, which are called non-global 
observables 14], the approximation of independent emissions, used in the analytical 
resummations, is not valid to single (which for some of these observables means leading) 
logarithmic accuracy. Non-global observables typically involve measurements of soft emis- 
sions over a limited part of phase space, a good example being energy flow distributions in 
a fixed rapidity-azimuth (rj — (ft) region. In fact in the case of the energy flow away from 
hard jets the function g\ in Eq. ( |1 . 1|) is absent (there being no collinear enhancement in 
the away-from-jet region). The leading logarithms in this case are thus single logarithms 
that are resummed in a function equivalent to 52 but this function cannot be completely 
calculated within an independent emission formalism. This is the case because the inde- 
pendent emission approximation of the QCD multi-parton emission pattern is strictly valid 
and intended for use in regions where successive emissions are strongly ordered in angle. 
The leading partonic configurations (those that give rise to the leading single-logarithms) 
for the away-from-jet energy flow are however those which include the region of emission 
angles of the same order in the parton cascade. Thus relevant single-logarithms also arise 
from multi-soft correlated emission which has been computed only numerically and in the 
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large- A^ c limit thus far fO| ||]. 

Since one of the main approximations used in analytical resummations, that of inde- 
pendent emission, has been shown to be inaccurate even to leading-logarithmic accuracy 
for some non- global observables like interjet energy flow, one is led to wonder about the 
leading-logarithmic accuracy that is claimed for parton showers in Monte Carlo event gen- 
erators, in these instances. The parton shower in HERWIG for instance relies on an 
evolution variable k 2 which in the soft limit is equivalent to ordering in angle [11, [Ef3]. 



Angular ordering of a soft partonic cascade, initiated by a hard leg, is a perfectly good 
approximation for azimuthally averaged quantities such as some e + e~ event shapes and in 
fact can be further reduced in these instances to an independent emission pattern, up to 
next-to-leading logarithmic accuracy. However, when looking at energy flow into limited 
angular intervals, one is no longer free to average soft emissions over the full range of angles, 
which means that one no longer obtains angular ordering at single-logarithmic accuracy. 
Thus one expects at least formally that the parton shower in HERWIG is not sufficient 
even to leading logarithmic accuracy for variables such as energy flow in inter-jet regions. 
The same statement should apply to the PYTHIA shower and even more strongly to ver- 
sions before 6.3 where the ordering variable is always taken as the virtuality or invariant 
mass and angular ordering imposed thereafter ||, which leads to insufficient phase-space 
for soft emission. Version 6.3 @ offers as an alternative the possibility to order the shower 
according to the transverse momentum of the radiated parton with respect to the emitter's 
direction (see || for more discussion on the transverse momentum definition), which yields 
a better implementation of angular ordering Q. We would like to point out that the ARI- 
ADNE Monte Carlo generator [17] has the correct large-angle soft gluon evolution pattern, 
which generates the non-global single logarithms in the large- N c limit. Since however the 
most commonly used and popular programs are the ones we mentioned before, we shall be 
interested in comparisons to the showers therein. 

This issue assumes some importance while considering for instance the tuning of the 
shower and non-perturbative parameters in Monte Carlo generators. If the tuning is per- 
formed by using data on a non-global observable such as energy flow away from jets one 
must at least be aware of what the accuracy is of the shower produced by the event genera- 
tor. If the accuracy is not even leading- logarithmic then one runs the risk of incorporating 
missing leading-logarithmic effects via tuned parameters. This situation is not optimal 
since, as far as possible, one would like to account only for subleading effects and incal- 
culable non-perturbative physics via the tuning. Moreover, the soft physics of non-global 
observables is not universal, the multi-soft correlated emission component being irrelevant 
in the case of global observables (those sensitive to soft emission over the full angular 
range). This difference in sensitivity to soft gluons, for different observables, would not be 
accounted for in case the non-global effects are tuned in once and for all. 

In the present paper we aim to investigate the numerical extent of the problem and 
to what extent non-global logarithms may be simulated by angular ordering and hence by 
parton shower Monte Carlo generators. In the following section we shall compare a fixed 
order 0(a 2 ) calculation of the leading non-global effect for energy flow into a rapidity slice 
with that from a model of the matrix element where we impose angular ordering. We shall 
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comment on the results obtained and in the following section examine what happens at all 
orders and whether our fixed-order observations can be extrapolated. Having compared 
the full non-global logarithmic resummation with its angular-ordered counterpart we then 
proceed to examine if our conclusions are borne out in actual Monte Carlo simulations. 
Thus we compare the results of resummation with those obtained from HERWIG and 
PYTHIA at parton level. This helps us arrive at our conclusions on the role of non-global 
effects while comparing Monte Carlo predictions to data on observables such as the energy 
flow between jets, which we report in the final section. 



2. Non-global logarithms vs angular ordering at leading order 

In order to explore the issues we have raised in the introduction, we pick the interjet energy 
flow (more precisely transverse energy Et flow) observable. Here there are no collinear 
singularities and the problem reduces to one where the leading logarithms encountered in 
the perturbative prediction are single-logarithms. While the nature of the hard-process 
is fairly immaterial in the large- N c limit to which we confine our discussions, it proves 
simplest to choose e + e~ — >2 jets and examine the Et flow in a chosen angular region. 
Given a phase-space region O, the Et flow is defined as 

E t = J2 E ti> (2- 1 ) 

where the sum runs over all hadrons (partons for our calculational purposes) and the 
observable we wish to study is 

1 r Qn da 

£(Q,Qn) = ~/ dE t —. (2.2) 
a Jo atj t 

The theoretical result for the integrated quantity £ was correctly computed to single- 



logarithmic accuracy in Ref. [14| and assumes the form 

S(Q, Qa) = exp{-4C F A n t] S(t), (2.3) 

where one has defined t 

[Q dk t a s (kt) , . Q f0M 

*( L )=/ T — o^~> L = ln 7T- ( 2 - 4 ) 

JQe~ L Kt 27T Qq 

The first factor in Eq. ( |2.3| ) above is essentially a Sudakov type term where Aq = J d-q^ 
represents the area of the region Aq. Note the colour factor Cp from which it should be 
clear that this term is related to multiple independent emission off the hard primary qq 
pair and in fact is just the exponential of the single- gluon emission result. 

The second factor S(t) is the correlated gluon emission contribution which starts with 
a term that goes as CfCao? s ln 2 (Q/Qn). This can be calculated fully analytically while 
the full resummed single-logarithmic calculation for S(t) is carried out numerically in the 
large N c limit. Before we turn to the all-orders result we aim to compare the analytical 
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leading-order computation with a model of the matrix element based on angular ordering. 
This will give us some insight into the issue at hand. 

In order to do so we start with the full matrix-element squared for energy-ordered two 
gluon emission from a qq dipole ab: 



M 2 (k u k 2 ) = AC F { " h) 



(aki)(bk 



C A (ah) | C A | C a \ (ab) 



2 (ak 2 )(k 1 k 2 ) 2 (bk 2 )(hk 2 ) \ 2 J (ak 2 )(bh 



(2.5) 



with the conventional notation (ab) = a-b with a, b and ki being the particle four-momenta. 
We define these four-vectors as below: 

0=1(1,0,0,1), (2.6) 

b = | (1,0, 0,-1) , 
ki = ktA (cosh 771 , cos (pi, sin <j>\, sinhryi) , 
k 2 = k t>2 (cosh?7 2 ,cos(/>2,sin(/>2,sinh77 2 ) , 

where Q is the centre-of-mass energy. 

We also separate the "independent emission" piece of the squared matrix element, 
proportional to C F , from the correlated emission piece proportional to C f Ca- 

M 2 (ki, k 2 ) = C 2 F W(k l )W(k 2 ) + C F C A W(h, k 2 ) . (2.7) 

It is this latter piece that is termed the non-global contribution at this order. 

We now wish to distinguish between a full calculation of the non-global contribution 
at O (ttj) and that based on an angular-ordered model of the squared matrix element. We 
first revisit the full result without angular ordering. Since only the CfCa piece of the 
result will be different in the angular-ordered approximation, we shall focus on this term. 
Using the momenta defined in Eq. ( [2.6D we obtain 

cosh (771 - T] 2 ) 



c F c A w(k u k 2 ) = 1 -^^ 



cosh (rji — 772) — cos( 



where we introduced the transverse-momentum fractions Xj = 2kfi/Q, and assume that 
x\ 3> x 2 , i.e. strong ordering of the transverse momenta. 

The non-global contribution is given by integrating the above result over the directions 
of the two gluons such that the softer gluon (^2) is in f2 while the harder gluon (k\) is 
outside, and over the scaled transverse momenta x% and x 2 . The integral over directions 
(including a phase space factor Q 4 /16) is given by 



Q 



CfC a ^ I d m ^ I dr/2-^ W(h,k 2 ). (2.9) 

Integrating over the energy fractions x\ and x 2 produces at the leading single-logarithmic 
level a factor -(1/2) \n 2 (Q/Q n )- The coefficient of the (f^) 2 \n 2 (Q/Q n ) term has a C F C A 
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or non-global contribution which reads 

cosh (rji - r/ 2 ) 



S 2 = -4C F C A I cftfcfi/ d m ^ 



cosh (rji — r/ 2 ) — cos( 



1 



(2.10) 



We now choose H as a slice in rapidity of width A77 which one can centre on rj = 
with its edges at rapidities —Arj/2 and Arj/2. We are free to take <f>i = and integrating 
over 02 gives the result 



dm / *?2 [C0th(772 - 77l)-l] , (2.11) 

-00 J~Ar)/2 



where we doubled the result of assuming 771 < 772 to account for the region rji > r/2- 

Now one is left with the integral over the gluon rapidities. In order to examine the 



main features of the final result, which were already elaborated in Ref. [13|, we introduce 
the rapidity difference y = f]2 — ill in terms of which one can reduce the above integral to 

J dyy (cothy - 1) + J dy At? (coth y - 1)J . (2.12) 

Let us concentrate on the case of a large slice where the result has an interesting behaviour. 
As one increases Arj the second integral in the sum above, from Arj to infinity, starts to 
become progressively less significant. This is because the integrand cothy — 1 rapidly 
approaches zero as y becomes large. The first term in the parentheses, on the other hand, 
gets its main contribution from the small y region. Its value as Arj — » 00 tends to 7r 2 /12. 
Thus what one observes as one increases Arj is that the contribution to the integral from 
Arj — > 00 starts to be negligible while the contribution of the integral from zero to Arj starts 
to become insensitive to its upper limit and hence the slice width A77, being dominated 
by the contribution from the small y region. This leads to a rapid saturation of the result 
as one increases Arj and the result quickly approaches 7r 2 /12. For instance the value at 
Arj = 2.5 is 0.818 while vr 2 /12 = 0.822. 

Now we recompute the above integral using an angular-ordered approximation of the 
squared matrix element. We expect that the angular ordering we introduce here should 
correspond to the contribution to the non-global logarithms that ought to be contained in 
Monte Carlo event generators based on angular ordering. The angular-ordered approxima- 
tion to the matrix element squared Eq. (|2.5| ) is obtained by modifying each dipole emission 
term therein as below: 

(ab) 1 — cos 6 a b 



(ak)(bk) lo 2 (1 — cos 9 a k) (1 — cos 9kb) 

_^ J_ / 9 (cosflqfc - cosflqfc) + 9 (cos 6kb - cos 9 ab ) \ 
to 2 \ 1 - cos 9 ak 1 - cos 9 k b ) 

where u refers to the energy of k. The second line above is actually equivalent to the full 
result if one can integrate freely over the azimuthal angles defined with respect to each of 
the legs of the emitting dipole, leaving a dependence on just the polar angles 9. However, 
since one places geometrical restrictions on the emissions k\ and k2, and in that respect 
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k\ has to be outside the gap while &2 inside, the azimuthal integration does not extend 
from zero to 2ir. The limits instead depend on the precise gap geometry. Ignoring this we 
wish to model the full matrix element squared by the angular pattern introduced above, 
corresponding to emission of soft gluons in well-defined cones around each hard emitting 
leg. 




An 

Figure 1: The coefficient of the leading order non-global contribution —S^ICfCa plotted as a 
function of the rapidity slice A?7 as given by both the full calculation and the angular-ordered 
approximation. The significant feature of saturation of the result for large slice-widths is apparent 
in both results. 



We note once more that the C% independent-emission term of the squared matrix 
element is left intact since the angular-ordered and full results are identical for this piece, 
as one would expect. Making the modification described in Eq. ( p. 13 ) in each term of the 
CfCa piece of the squared matrix element Eq. (^5|) and integrating over gluon directions 
we obtain the coefficient 52 in the "angular ordered" (AO) approximation. We plot the 
numerical result in this approximation as a function of the gap size in Fig. |l| along with the 
full result. One can immediately observe that for small gap sizes the AO and full results 
are essentially identical. As one increases the gap size one notes a numerically significant 
difference between the two results although this is at best moderate. For instance for a 
slice of width Ar] = 2.5 one observes that the AO result is lower by 10.67% than the full 
result. Additionally it is interesting to observe that the notable feature of saturation of S2 
for a large gap size is preserved by the AO approximation. 

The reason the saturation property is preserved is because, as explained previously 
in detail, it arises from the region where the two gluons (respectively in and outside the 
gap region) are close in angle or equivalently from the region of integration r]\ — r}2 < Ary. 
Moreover, the bulk of the non-global contribution for any gap size arises from the region 
where the emission angles of the two soft gluons are of the same order. The contribution 
from configurations with the softest gluon at large angle relative to the next-softest gluon 
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are small and vanish rapidly as we make the rapidity separation 771 — 772 large. 

In the AO approximation one requires the softest gluon &2 to be emitted in a cone 
around the hard emitters k\ and either the emitting quark leg a or b, depending on whether 
one is looking at emission by dipole ak\ or bk\. The size of the cone is equal to the 
dipole opening angle. Thus the important region where k\ and ki are collinear is perfectly 
described by the AO model. Only the region where k<i is emitted at an angle larger than the 
cone opening angle would not be covered in the AO approximation and the contribution of 
such a region should be relatively small as we observe numerically. We mention in passing 
that these conclusions described explicitly for a rapidity slice are expected to hold for a 
general gap geometry and we explicitly checked the case of a square patch A77 = A<f> in 
rapidity and azimuth. 

In the following section we shall examine the impact of the AO approximation at 
all orders to determine whether the encouraging fixed-order finding, that an AO model 
reproduces the characteristics and is numerically reasonably close to the full non-global 
result, can be extended to all orders, as one may now expect. 



3. AO approximation at all orders 

We now study the AO approximation by using the large N c evolution algorithm that was 



described in Ref. [ 13 ] , suitably modifying it to take account of the angular ordering re- 
quirement. This should enable us to estimate how non-global logarithms will be simulated 
in an angular-ordered parton shower event generator. The algorithm works as follows. 
To compute the non-global contribution S(a s L) where L = ln(Q/Qn) one considers the 
probability Pc(L>) of a configuration C that does not resolve gluons above scale L, in other 
words those with energies below Qe~ L . The evolution of this configuration to another con- 
figuration C with larger resolution scale V or equivalently smaller energy scale, proceeds 
via soft emission of an extra gluon k! from the configuration C: 

P c ,(L') = a s (L')A c (L,L')Pc(L)F c (8 f ,cj ) '), (3.1) 

where Ac(L, L') represents the summation of only virtual gluons between the scales L and 
L' , Fc(0',4>') represents the angular pattern of emission of gluon k' from the system of 
dipoles in the configuration C and a s = a s /(2ir). One has explicitly 

Fc{dk,9k)= y. jz 7. r-jz 5— r. (3.2) 

^ .. (1 - cos e ik ) (1 - cos e jk ) 

dipoles— ij J 

The same dipole angular pattern enters the pure virtual evolution probability (or form 
factor): 

lnA c (L,L') = -J dL" j d cos 6 d(f>a s (L")Fc (6, 0). (3.3) 

The probability that the interjet region Q stays free of real emissions below a given 
scale L, is then given by summing over corresponding dipole configurations: 

£(Q,Qn)= Yl P ^ L )- ( 3 - 4 ) 

C I O empty 
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In order to obtain our angular-ordered results we need to modify the angular emission 
pattern Fq, as before for the fixed-order case, so we define 



Fc(O k , A) ao ^ £ 2cJ Q (™ e *-™ e «K ei ™ 9ik -™ eii) ). (3.5) 

,. , .. V l-COS0jfc l-COSfc'jfc J 

dipoles— ij 

Making the replacement Fc(0k, ^fe) — * Fc{@ki <t>k)AO ° ne modifies both real and virtual 
terms and obtains the result from our angular-ordered model at all orders: 

Z AO (Q,Qn)= £ p c,A0(L). (3.6) 

C | f2 empty 

Having obtained Sao we can compare it with the full result without angular ordering. 
In Fig. |2] we plot the full and AO results for X(£) as a function of t, for a slice of unit width 
A77 = 1. 




Figure 2: The integrated cross-section as a function of t in the full calculation and the AO 
approximation. The primary result is also shown for reference. 



One notes the relatively minor difference between the full and the AO results which 
indicates that the contribution to the full answer from regions where one can employ 
angular ordering, is the dominant contribution. For the sake of illustration we focus on the 
value t = 0.15 which corresponds to a soft scale Qn = 1.0 GeV for a hard scale Q = 100 
GeV. For the rapidity slice of unit width we note that the result for Sao(*) is 9-68 % higher 
than the full result. At the same value of t, the difference between the full and the primary 
result, i.e. exp(— 4,Cpt), is around 75 %, thus indicating that the AO approximation is 
much less significant than the role of the non-global component itself. Similar observations 
hold regardless of slice width. 
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One can also directly study the impact of the AO approximation on the pure non-global 
contribution S(t). The primary contribution is unaffected by angular-ordering and can be 
divided out from the result for Eao(^) to give us SaoW- We first take the example where U 
is a rapidity slice and consider different values for the slice width Ar]. We illustrate in Fig. || 
three choices for the slice width Arj = 1.0, 2.0, 3.0 with the full non-global contribution 
S(t) and that in the AO model. We note that in both full and AO cases the feature of rough 
independence on the slice width Ar] is seen, as one can expect for sufficiently large slices. 
The AO curves are somewhat higher than the full ones indicating a somewhat smaller 
suppression than that yielded by the full calculation. 




0.05 0.1 0.15 0.2 0.25 



t 

Figure 3: The resummed non-global contribution S(t) as a function of t in the full calculation 
and the AO approximation for different values of the slice width Ar]. The upper set of curves 
correspond to the AO case and reflect that in that approximation a slightly smaller suppression is 
obtained than from the full calculation corresponding to the lower set of curves. The feature of 
rough independence on the slice width Ar] is visible in the full case and is preserved by the AO 
approximation. 



Similar studies can be carried out for different geometries of f2. For a square patch in 
rapidity and azimuth with Ar] = A(f) = 2.0, the full and angular-ordered results for 
are shown in Fig. ||. The difference is seen to be small over a wide range of t. Once again 
focusing on the t = 0.15 value, one notes that the AO approximation is only about three 
percent above the full result. At t = 0.2 this difference rises to 9.75 %. Corresponding 
results for S(t) for the same square patch, obtained by dividing by the primary result, are 
plotted in Fig. || and once again only a small to moderate effect is observed over the t range 
shown. 



We have thus observed that modifying the evolution code [ 13 1 used to compute the non- 
global logarithms, to impose angular ordering on them, only has a moderate effect on the 
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t 



Figure 4: E(£) vs t for a square patch in rapidity and azimuth, Aij = A<f> = 2.0. Primary, full and 
angular-ordered (AO) curves are shown. 




t 



Figure 5: The non-global contribution S(t) as a function of t for a square patch in rapidity and 
azimuth, A77 = Acf> = 2.0. 

quantity S(t). This effect becomes even less significant for the quantity £(t) = Sp(t) S(t) 
since the primary contribution Xp(i) = exp[— 4Ci?Ar>i] is unchanged by imposing angular 
ordering, which we also explicitly checked with the code. 



- 12 - 



Having thus noted the small effect of the AO approximation within our model we 
would not expect much difference, in principle, between the results from an event generator 
based on angular-ordering in the soft limit (HERWIG) and the full non-global results. 
For PYTHIA, prior to the version 6.3 one may expect to see differences since angular 
ordering was imposed on top of ordering in the virtuality (invariant mass) of a splitting 
parton which leads to known problems with soft-gluon distributions, as discussed in ||. 
In Ref. [18], where colour coherence effects were observed and studied at the Tevatron 
collider, it was in fact found that, unlike HERWIG, the PYTHIA event generator was not 
able to acceptably reproduce experimental observables sensitive to angular ordering. One 
may expect, however, that the new PYTHIA model |, § (where, to our understanding, the 
improved shower, ordered in transverse momentum, better accounts for angular-ordering) 
results comparable to those from HERWIG may be obtained. In the next section our aim 
is to explore these issues and see if our expectations, outlined above, are indeed borne out. 



4. Comparison with HERWIG and PYTHIA 

In this section we shall focus on actual comparisons to results from HERWIG and PYTHIA. 
In order to meaningfully compare the results of a leading-log resummation with the parton 
level MC results, it is necessary to minimise the impact on the MC results of formally 
subleading and non-perturbative effects that are beyond full control and hence spurious. 

In order to suppress subleading effects one needs to carry out the comparisons to the 
MC generators at extremely high values of Q, and hence we chose 10 5 GeV. Thus effects 
that are formally of relative order a s (Q) or higher can be expected to be negligible. A sign 
of this is the fact that at such large Q values the MC results one obtains do not depend 
on Q other than via the single logarithmic variable t for a large range of t. It is clear that 
such high Q values are beyond the reach of current or imminent collider experiments but 
since we are interested only in the dependence on t, the Q value is fairly immaterial for our 
purposes. In fact one can take the conclusions we make for a particular t value at Q = 10 5 
GeV and translate that into a value of Qn for an experimentally realistic value of Q. 

A clear source of uncertainty in this procedure is the different definitions of a s in the 
resummation and the MC programs. In all resummed predictions we have used the LL 
expression for t: 

,_ 1 , 1 R _ nC A - AT R n f 

'"44 l-2a s (Q)/|L' P0 ~ 12vr ' ^ 

with Po corresponding to nj = 6 and L given in Eq. ( |2.4[ ). The coupling a s (Q) is in the MS 
scheme, and is obtained via a two-loop evolution with 6 active flavours from the input value 
a s (Mz) = 0.118. This is to ensure that the resummed prediction is a function of a s (Q)L 
only. HERWIG instead exploits a two-loop coupling in the physical CMW scheme |12| | 
with a s {Mz) = 0.116, while PYTHIA uses a one-loop coupling corresponding to a s (Mz) = 
0.127 0| . The values of t corresponding to different definitions of a s (computed according 
to Eq. ( |2.4f ) ) are found to be compatible within 10% in the considered Et range. This does 
not lead to appreciable modifications in the resummed curves plotted as a function of Et 
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rather than t, in this section. Thus the comparisons we make below to the Monte Carlo 
results at a particular value of Et are not significantly affected by the issue of the somewhat 
different definitions employed in the resummation and the various Monte Carlo programs. 

Another effect, not accounted for in the resummation, that is potentially significant, 
is the effect of quark masses (which would arise due to excitation of all flavours). These 
effects however can be safely neglected at the value of Q we choose. In particular we also 
note that the presented MC curves are obtained by allowing the top quark to decay, but 
we have explicitly checked that we obtain almost identical results if we force the top quark 
to be stable. 

With the above observations in place, we start with the comparison to HERWIG which 
has a parton shower which is ordered (in the soft limit) in angle and thus one would expect 
results in line with those obtained via our AO model, introduced in previous sections. In 
Fig. H we show the results obtained from HERWIG compared to those from resummation for 
a rapidity interval of unit width. We note here that in order to obtain a sensible behaviour 
for the resummed predictions at large Et, it was necessary to match the resummed results 
to exact fixed-order estimates. We carried out the so-called log-i? matching |l| to both 
leading and next-to-leading order (obtained from the numerical program EVENT2 
but at the values of Et we have shown here, no significant difference was observed. The 
curves plotted in Fig. ^ are matched to NLO while the HERWIG results contain matrix- 
element corrections pcfl . We observe that a very good agreement between HERWIG and 
the full and AO curves is seen over a significant range of Et values. We have also included 
the value of the variable t as a function of Et to enable us to extrapolate our conclusions 
to lower centre-of-mass energies. 

The comparison to PYTHIA is shown in Fig. g. We use version 6.3 and consider the 
old model, with showers ordered in virtuality and forced angular ordering, as well as the 
new model, where the emissions are ordered in transverse momentum. We note that the 
results obtained from PYTHIA with the new parton shower appear to be in reasonable 
agreement with the resummed curves including non-global logarithms, the situation being 
comparable to the quality of agreement one obtains with HERWIG. The same is not true 
for the old PYTHIA shower and a significant disagreement between the result there and 
the resummed curves is clearly visible. 

In order to be more quantitative we focus on Et = 10 GeV which corresponds to a value 
of t = 0.15. Here we note that the difference from the full resummed curve is respectively 
for HERWIG, PYTHIA (new) and PYTHIA (old) approximately -10%, +7.5% and 50%. 
The difference between a resummed primary contribution and the full non-global result 
is, at the same value of Et, 25%. We would then infer that if a variable of this type is 
chosen to tune for instance PYTHIA with the old shower (with ordering in the invariant 
mass) one includes potentially as much as 50% of the leading-logarithmic perturbatively 
calculable contribution, to model-dependent parameters and incalculable effects such as 
hadronisation and the underlying event. 

We have carried out our study for slices of different widths and obtain comparisons 
with HERWIG that are generally satisfactory. The same appears to be true of the new 
PYTHIA algorithm but here problems seem to crop up as one increases the slice rapidity. 
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Figure 6: The distribution a 1 da/dE t for a slice of Arj = 1 and Q = 10 5 GeV compared to parton 
shower results from HERWIG. 

In Fig. |8] we present the comparison with both HERWIG and PYTHIA, but for a slice 
width Ar/ = 3.0. We observe that for a larger slice the new PYTHIA shower at lower Et 
values yields a result that is significantly below all other predictions. The reason for this is 
not entirely obvious to us and we would welcome further insight into this observation. We 
have also carried out studies at other intermediate slice widths e.g. Arj = 2.0 and it appears 
that the new PYTHIA curve starts to deviate from the resummed results at a value that 
is exponentially related to the slice width. This may signal that the new ordering variable 
in PYTHIA is perhaps not entirely satisfactory at large rapidities but as we mentioned a 
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Figure 7: The distribution a 1 da/dE t for a slice of A77 = 1 and Q = 10 5 GeV compared to parton 
shower results from PYTHIA. 

more detailed study is required to draw firm conclusions on this issue. 
5. Conclusions 

In this paper we have examined the role played by angular ordering in the calculation 
of the leading single-logarithmic terms that arise for non-global observables such as the 
away-from-jet energy flow. While it has been clear for some time [[0|, 14 1 that the fully 



correct single logarithmic resummed result cannot be obtained via use of angular ordering 
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Figure 8: The distribution a^dajdEt for a slice of A 77 = 3.0 and Q = 10 5 GeV. 



the question remained as to how much of the full result for such an observable may be 
captured by using the approximation of angular ordering. The reason this question arises in 
the first place is mainly because angular ordered parton showers are employed for example 
in Monte Carlo event generators such as HERWIG. Given the importance of these event 
generators as physics tools it is vital to understand the accuracy of the different ingredients 
thereof (such as the parton shower). 

While the accuracy of the parton showers is generally claimed to be at least leading- 
logarithmic, this statement ought to apply only to those observables where the leading 
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logarithms are double logarithms, i.e. both soft and collinear enhanced. However, to the 
best of our knowledge, there has been no discussion yet in the literature about non-global 
observables where the leading logarithms may be single logarithms instead of double log- 
arithms and the accuracy of the parton showers in such instances. Since observables of 
the type we discuss here (energy or particle flows in limited regions of phase space) are 
often used in order to tune the parameters of the Monte Carlo algorithms (see e.g. [21] for 
examples and references), it is important to be at least aware of the fact the perturbative 
description yielded by the parton shower, may in these cases be significantly poorer than 
that obtained for instance for global observables. We have thus chosen one such observ- 
able and carried out a detailed study both of the role of angular ordering as well as the 
description provided by the most commonly used Monte Carlo event generators HERWIG 
and PYTHIA, compared to the full single- logarithmic result (in the large N c limit). 

We find that in all the cases we studied, involving energy flow into rapidity slices or 
patches in rapidity and azimuth, angular ordering captures the bulk of the leading logarith- 
mic contribution. This is a comforting finding but there remains the issue of precisely how 
angular ordering is embedded in the parton shower evolution for HERWIG and PYTHIA. 

For HERWIG where the evolution variable in the soft limit is the emission angle one 
expects the agreement between parton shower and the leading-log resummed descriptions 
to be reasonable and we find that this is in fact the case. 

In the case of the PYTHIA shower (prior to version 6.3) angular ordering is imple- 
mented by rejecting non- angular-ordered configurations in a shower ordered in virtuality. 
In this case it is clear that the description of soft gluons at large angles will be inadequate 
H and this feature emerges in our studies. From this we note that a discrepancy of around 
50% could result while comparing PYTHIA to the correct leading-log result. This differ- 
ence would be accounted for while tuning the parameters of PYTHIA to data and must be 
borne in mind, for instance while making statements on the tuning of the hadronisation 
corrections and the underlying event into the PYTHIA model. This is because a tuning 
to energy flows would mean that significant leading-logarithmic (perturbatively calculable) 
physics is mixed with model-dependent non-perturbative effects which does not allow for 
the best possible description of either. Moreover, the non-global effects are not universal 
and thus incorporating them into the generic shower and non-perturbative parameters will 
lead to a potentially spurious description of other (global) observables. 

The new PYTHIA shower, ordered in transverse momentum and with a more accurate 
treatment of angular ordering, does however give a good description of the leading loga- 
rithmic perturbative physics, comparable to that obtained from HERWIG. However, for 
large rapidity slices we find that problems emerge in the description provided by PYTHIA 
even with the new shower. The origin of these problems is not entirely clear to us and we 
would welcome further insight here. Hence, we strongly emphasise the need to compare 
the shower results from HERWIG and PYTHIA while carrying out studies of observables 
that involve energy flow into limited regions of phase space. Where this difference is seen 
to be large, care must be taken about inferences drawn from these studies about the role 
of non-perturbative effects, such as hadronisation and the underlying event. We believe 
that further studies and discussions of the issues we have raised here are important in the 
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context of improving, or at the very least understanding, the accuracy of some aspects of 
Monte Carlo based physics studies. 
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